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Abstract 

We study through symmetry principles the form of the functions in the generalized scalar-tensor theories under the self- 
similar hypothesis. The results obtained are absolutely general and valid for all the Bianchi models and the flat FRW one. We 
study the concrete example of the Kantowsky-Sach model finding some new exact self-similar solutions. 

1 Introduction 

Current observations of the large scale cosmic microwave background suggest to us that our physical universe is expanding 
in an accelerated way, isotropic and homogeneous models with a positive cosmological constant. The analysis of Cosmic Mi- 
crowave Background (CMB) fluctuations could confirm this picture. But other analyses reveal some inconsistencies. Analysis 
of WMAP data sets shows us that the universe might have a preferred direction. For this reason, it may be interesting to study 
Bianchi models since these models may describe such anisotropies. 

The observed location of the first acoustic peak of the temperature fluctuations on the CMB corroborated by the data 
obtained in different experiments [1], indicates that the universe is dominated by an unidentified "dark energy" and suggests 
that this unidentified dark energy has a negative pressure [2]. This last characteristic of the dark energy points to the vacuum 
energy or cosmological constant as a possible candidate for dark energy. Although, it is a general belief that the current 
curvature of the universe is negligible and mostly the universe is considered with a flat geometry, recent observations support 
the possibility of a non-flat universe and detect a small deviation from k = [3]. For example, evidence from CMB and also 
supernova measurements of the cubic correction to the luminosity distance favour a positively curved universe [4, 5]. 

In order to explain the current acceleration of the universe, within General Relativity (GR), it is necessary to introduce a 
new type of energy with a negative pressure. Between the different possible approaches is one which consists of considering 
so-called dark energy (DE). There are several candidates for DE, where the simplest one is the cosmological term A. However, 
this choice has several drawbacks such as coincidence and fine tuning problems. For this reason other models have been 
proposed. Examples of such models are quintessence [6, 7], K-essence [8], or chameleonic fields in which a scalar field is 
coupled to matter [9], etc. Hence it is natural and important to consider a variable A— term in more general frameworks 
where, furthermore, other quantities, such as the Newton gravitational constant may be considered as dynamical. One such 
class of theories are the scalar-tensor theories (STT) of gravity. This class of models has received a renewed interest in recent 
times, for two main reasons: First, the new inflationary scenario as the extended inflation has a scalar field that solves several 
problems present in the old theories. Secondly, string theories and other unified theories contain a scalar field which plays 
a similar role to the scalar field of the STT. The scalar-tensor theories started with the work of R Jordan in 1950 [10]. A 
prototype of such models was proposed by Brans and Dicke in 1 96 1 [11]. Their aim for presenting this model was to modify 
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Einstein's theory in such a way as to incorporate the so called "Mach's principle". These theories have been generalized by 
P.G. Bergmann [12], K. Nordtverdt [13] and R. T. Wagoner [14]. For a recent review of this class of theories we refer to [15] 
and [16]. 

In this paper we want to consider a family of scalar-tensor theories with a dynamical cosmological constant [17] and with 
a potential [15], that is equivalent to a time dependent cosmological constant. Recently several authors have considered the 
cosmological consequences of a time varying cosmological constant. Most of them introduce the time dependence in an ad 
hoc manner. In this work we consider an equivalent problem in the well known general scalar-tensor theory of gravity where 
the time dependence can occur in a natural way, without any new assumption or modification of the theory and provide an 
explanation for the acceleration of the universe expansion [18]. 

Many authors have studied these general scalar-tensor theories. They use the observational data in order to obtain restric- 
tions or constraint between the functions that appear in the action to obtain an accelerated model [19] -[22]. Our approach is 
different: We want to derive these functions from symmetry principles as self-similarity. We shall carry out our study under 
this assumption and state some general theorems that are valid for all the Bianchi models and of course for the flat FRW one. It 
is most appropriate for us to work in the Jordan frame (JF), in which the physical quantities are those that are being measured 
in experiments, even though the Einstein frame (EF) often provides a better mathematical insight. 

The study of self-similar (SS) models is quite important since a large class of orthogonal spatially homogeneous models 
are asymptotically self-similar at the initial singularity and are approximated by exact perfect fluid or vacuum self-similar 
power-law models. Exact self-similar power-law models can also approximate general Bianchi models at intermediate stages 
of their evolution. This last point is of particular importance in relating Bianchi models to the real Universe. At the same time, 
self-similar solutions can describe the behaviour of Bianchi models at late times i.e. as t — > °° [23]. 

This paper is organized as follows. In section two we start by considering a particular formulation of the theory. In this 
case the cosmological constant is introduced directly by the function A(0) [17]. We study this model through two different 
approaches. The first consists of studying the effective stress-energy tensor under the matter collineation approach. This 
method allows us to obtain relationships between the physical quantities as well as to determine the exact form of the scalar 
field . The second approach consists of studying the wave equation under the Lie group method. By imposing a particular 
symmetry we are able to determine the exact form for each of the functions that appear in this equation. We summarize 
all the results by stating a very general theorem. In section three we study a very general scalar-tensor theory. In this case 
the cosmological constant is introduced by the potential. Following the same exposed procedure as in the above section, we 
are able to determine the exact form that all of the unknowns involved in this model must have. We show, from the stated 
theorem, how different versions of this theory arise, which are the standard Brans-Dicke theory, the induced gravity model 
[24] and a very particular solution where the effective gravitational function is constant. In section four we study a chameleon 
Jordan-Brans-Dicke model. In order to show how all the obtained results work, in section five we study a particular example 
which is the Kantowski-Sach model. We start this section by showing that this metric admits a homothetic vector fields and 
then we study several models. We put special emphasis on comparing the solutions obtained in each case. In section six we 
end by summarizing all the results. We have added an appendix where we study in detail one of the equations obtained in 
section 2 in order to show with out any doubt that the Brans-Dicke parameter co (<p) must be constant in this framework of 
self-similar solutions. 

2 Cosmological models with dynamical A in scalar-tensor theories 

Following to Will (see [17]) we start with the action for the most general scalar-tensor theory of gravitation 



S = 




g <j)R 



<0<0:g' i O J j 
4> 



+ 2</>A(0) +S NG , 



(1) 
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where g — det(gy), G* is Newton's constant, Sng is the action for the nongravitational matter. We use the signature 
( — ,+,+,+). The arbitrary functions 0)(<j>) and A(0) distinguish the different scalar-tensor theories of gravitation, A(0) 
is a potential function and plays the role of a cosmological constant, and (0(<j)) is the coupling function of the particular 
theory. 

The explicit field equations are 

fly-igy*=J^7y+A(0)gy + ^ , (2) 

dco , , i _ , / , dA 



(3 + 2co(0))D0 =8tzT- — f ,f -20 U— -A(0)J , (3) 
where T = T? is the trace of the stress-energy tensor. The gravitational coupling G s g(t) is given by 

„ , . /2co + 4\ G* 



2.1 Matter collineations 

We may calculate the relationship between the quantities (in a SS approach) by calculating the matter collineations. There- 
fore we have to compute (we use unit where 8tt = c = 1 ) 



',/f " ^uf /f ') + £ (0 ;u Stf O ) +A((j>) 8i j, (5) 
' ( 7}f ) 



neff_ 1 T , »0) 
f 



where //G stands for a homothetic vector field. For simplicity we have used a fiat FRW metric but we would like to emphasize 
that all the obtained results are absolutely valid for all the Bianchi models, since we only look for the behaviour of the physical 
quantities instead of restriction on the scale factors. Therefore the homothetic vector field (HFV) is 

HO = (t + t Q ) d, + ( 1 - (r + 1 ) H) xd x + ( 1 - (r + 1 ) H) yd y + ( 1 - (t + r ) H) zd z . (7) 

Note the non-singular character of the HVF, nevertheless for simplicity in the calculations we use the singular case. 

1. Ti=^-%, 

^Ho(^Tij)=0, -fp0' + rp'0 + 2p0 =0, (8) 





obtaining 

p' 0' _ 2 p 

P t 



= r 2 . (9) 



2. t 2 = (o($)4>- 2 hat, & l 

^°(^(^"5 ft ^'')) =l) ^ ttf(a^- 2(0) + 2t co0„0 +20)0,0=0, (10) 
and therefore 

2 ^ «00 ^ 0, 

or 
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^0,-2^ + 2^ = -- ^_ 2 ^ + 2 i, = _2 O)(0)%=r- (12) 

co n r r co f r vv > 2 
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note that ft)' = (O^tj),. 

We find the next solution for Eq. (1 1) 



d<p-C x lnf + C 2 =0. 



(13) 



For example 

• co (cp) = conf 



/ — dcp — klnij), 
J (p 



= W. 



(14) 



this is the unique solution mathematically possible (compatible with the SS hypothesis, see the appendix for an 
explanation). 

• <o(p) = <p° 

r<t> <W 2 2 „,, 2 



/ !L- d(p = -r /2 

J (p a 



^ a / 2 = dln? 0=Ci(lnf) 2/ ". 



If ft) = const, then 



f0 f 2 -/0 ft 0-0 t = 0, 



and therefore we get the following ODE 



? 0, 



=exp(-C 3 -C 2 lnf) = o f C2 



3. r 3 = 0- 1 (^-g y n« 



i.e. 



2Ho[-j{<l>y-8iPt) ) =0, 



f («>„/ + 0,/" - jf - 0,y ) + 20,/ 



2// - — ^ 0" - 2//— + 2H V 
/ 



2(f'-2fH) = 0, 



(15) 
(16) 

(17) 

(18) 

(19) 
(20) 



note that H = ht 1 , /1 £ R + . These equation are different for each Bianchi model and we only obtain restriction on the 
scale factors. 



4. T 4 =A(<j>) gij , 



i.e. 



fA^0'+2A = O 



if HO (A(0)gy)=O 



A, 2 
A V f 



A = A r 



-2 



(21) 
(22) 



where A' = A^0'. 



2.2 Lie groups 

We are going to study the Eq. (3) through the LG method, i.e. we study the kind of functions A(0) and ft) (0) such that 
this equation is integrable. We start by rewriting it in an appropriate way 



(3 +2ft)(0)) (0"+fe-y) = Cr a +B (A - 0A^) - 0' 2 ft),!,, 



(23) 
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where /; = const., h S R + , B = 2, and C = 87r(l — 3y)p . Note that we are taking into account the conservation equation 
divT = 0, i.e. p = p t~ a , where a = (1 + y)/i, and// = ht~~ l . 

We need to solve the following system of PDE 

mfa-W$n = 0, (24) 
2/ 1 r- l W^+W77 # -2W^,-(2W- l ^-a) W )7 7 + ^7 7 ^ = 0, (25) 

-3 (2ty (A - 0A,) +Cf-«) ^ +/ 1f - 2 W(^ r - 1) +2W77 (0 - 

-W% tt + 2hr 1 (0 <t , (1 -(3 + 2co)W" 1 )77+2co (| ,77 ; =0, (26) 
[B (0 2 A W + - A) + 2aty W _1 (Bp (A - 0A^ ) + Cr a )] rj - 2 [Bp (A - 0A* ) - Cr a ] ^ + 

+ aCr 0_1 § + [50 (A - pA^,) +Cr a ] ri^+ht^Writ+Wr),, = 0, (27) 
where W = (3 + 2co (</>)) . notice that (2fi) + 3)W -1 = 1, so Eq. (26) yields 

- 3 (Bp (A - 0A ) + Cr a ) ^ + / 2 r 2 W (f |, - §) + 2Wri t ^ - W^ tt + 20*17, = 0. (28) 

The symmetry % =t,7] = np, brings us to obtain the following restriction on A (p) . From Eq. (27) we get 
[B (p 2 A H + pA* -A) +2co^W- 1 (Bp (A - pA^,) + Cr a )] np - 2 [Bp (A - pA$) - Cr a ] + 

+ aCr a + [Bp (A -pA^)+ Cr a ] n = 0, (29) 

so 

Cr a (a+n-2 + 2n(i)^pW~ l ) = 0, (30) 

and 

n (p 2 A H + pA* - A) + (A - pA$) (-a) = 0, (31) 
where we have taken into account Eq. (30) therefore we have obtained the following ODE for A (p) 

'n + a\ ( Aa A 



</> P 2 



and whose general solution is 



(32) 



A(p)=AoP^" + ^ + -^-p, (33) 
Zn + a 

we choose 

A(p) =A <^ ( " +a) = A f- ( " +a) , (34) 
since from our result from the MC approach we already know that, A (f ) = t ~ 2 , so this means that n + a = 2. 

In the same way we may calculate the restriction on co (p) , 

, 2-a-n . (2-a-n\ (3 + 2co) 

whose solutions are 

co(p) = pir(n+ a - 2 ) e - 2 " c i-^ = (o n t- {n+a - 2) -^, »„=const, (36) 

therefore we obtain 

(0 (p) = const., n + a = 2. (37) 

In an alternative way, from Eq. (25) we get 

2-a-2n\ CO* m f^-r>. , r>.*h(2-n-a) _ m JZ-n-a 



«>W = ( Jy, => co(P)=C & +C 9 p^ z -"- a > = co t z -"- a , (38) 

and therefore 

co(p) = ©o = const, (39) 
since 2 = « + a. In the appendix we shall give and alternative and detailed proof of this result. 
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Theorem 1 The scaling and in particular the self-similar solution admitted for the FE (2-3) have the following form 

= *o(' + 'o)", A WO = A <^ ( " +a) = A Q (t + t Q )-^ +a) , 
with n + a = 2, therefore A (?) = Ao (f + ?o)~ 2 ■ The Brans-Dicke parameter is constant 

co (<p) = const. 

fl«(ip = p (t + t y a , a = (l + Y)h. 

3 The General case 

We start by defining the action [21] 



and therefore the FE read 



FG llv = SnT^ v a " er +Z 



[F-^F-y - g^UF] - U[(j))g^ 



and 



2ZU$ =F$R-Z^ a -2U$ , 



(40) 



(41) 



(42) 



where R is the scalar curvature. 



The effective gravitational constant G e g between two test masses measured in laboratory Cavendish-type experiments is 
given by 



Geff : 



2Z(<j>)F+4Fi 
2Z(0)F + 3F 2 



(43) 



3.1 Matter collienations 



We may define 



and therefore, as above: 



(44) 



1. 7i =F(0) 



so 



-%0 FWTy =0 



F Y p t 



2. r 2 = ^(^-^ y ^) 



HO 



F(</>) 



thus 



Z y F r 



tpF^<j)'-tp'F-2pF = 0, 



P -2 

F =t ' 



-f 2 = r 2 . 



(45) 
(46) 



tpFZ^ a - tZF^' 2 + 2fZF</>" + 2ZF0' = 0, (47) 



(48) 
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^' ^ 3 = 7W) (^M^; v — Suv^F) depends on the metric so we only obtain restriction on the scale factors. 
AT- a 

^ho (jF^ySfvJ = t<j)%F-tUF^' + 2<j)U = 0, (49) 

and therefore 

= <=► " (50) 



3.2 Lie groups 

For this model we have to solve the following equation 

2ZD(j)=F^R-Z^ a -2U < p, (51) 

that we may rewrite in the following form 

2Z(0 ff +to~y) =Cr%-Z^ a -2U^ (52) 
where we have assumed <j) = (/), and the derivatives respect f are denoted by a comma. Note that R w Cf ~ 2 , with CsR. 
The standard procedure brings us to outline the following system of PDE: 

(ZZ H -Zj) 77 + ZZ^ +4Z 2 hr 1 ^ +2Z 2 77 W -4Z 2 <^ = 0, (54) 
3(2V^-Cr 2 F 1 ,)^+2hr 2 Z(t^-^)+4Zr lt4l -2Z^ l +2Z^r lt = 0, (55) 

[Cr 2 (Zfy-ZF^) +2(ZU„-ZM)] ?7 +Z(Cr%-2U^) 77^ + 

2Z(2f/ ()) -Cr 2 />)^+2Cf- 3 Z J F^+2Z 2 (/!f- 1 77, + 77„) =0. (56) 

The symmetry ^ = t , 77 = «0, brings us to obtain the following restrictions. From Eq. (54) 

Z 2 77 Z 2 Z 
Zm = ~2~Zq — > ==> = — — , => Z(<j>) = Zo<j)~"\ (57) 

where ?7i € R. From Eq. (56) we get 

2(ZU^-Z^)r]-2ZU^^+AZU^ t = 0, (58) 
Cr 2 [(Z^-ZF w ) 77 + ZF^-2ZF^ + 2r~ 1 Z/>4] = 0, (59) 



and therefore 



so 



while 



obtaining 



Za fn-2\U (n-2 \ U* 

u *i + \—)i> = _w )f' (60) 

f/(^) = C2 + f/o^ _ " ( " m_2 ' ,+2) , (61) 

F W= F *y + y> F w -(l-m)^, (62) 

F(0)=C 1 +F O (/. 2 -'". (63) 
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Theorem 2 The scaling and in particular the self-similar solution admitted for the FE (41-42) have the following form 

= O (r + r )" , Z(4>)=Z Q r m = Z Q (t + t Q )-'" n , (64) 

with 



while 



p=p (t + to) °\ a = (l + y)h, 
F(<j))=C l +F (j) 2 - m , U (<j))= C 2 + U $->"- 2n+ V 



(65) 
(66) 



This theorem state that we have a set of theories which admit scaling and in particular self-similar solutions if the involved 
functions take this particular form. The action reads 



s= hl dAx ^{\ [^ m *-r m <M a -2E/w] +^m| 



(67) 



So setting different values for the constant m we obtain different theories. For example if m = (induced gravity case) then 
the action (40) yields 



0=0 O (? + * O )'\ Z(#)=Zo, F(<p) = F^\ C/(0) = t/ o 0" (2 ^ 2) , 
and if m = 1 (usual JBD theory with a potential) then the action (40) yields 



G e ff « 0" 



i - 

%TiJ V S \2 



ft) 



J R--0, Cf a -2f/(0) 



= o (f +?<))*, z((j>)=z Q ^ 



l.e. 



Z(0) 



F(0)=F o 
<»(0) 

' 



t/(0) = t/ o 0" ( "~ 2) , 
ft) (0) = const. 



G e ff « 0~ 



and to end if m = 2, then the action (40) yields 



1 



S = — / <f ji\ 
8?r , 



to 



0" 



we get: 



= o (f + f o )", Z(0)=Z o 0- 2 , F(0)=F o , U(<j>) = U <j>^, Geff « const. 



This particular case is very similar to the scalar field cosmological model with G e g « const. Notice that this model is quite 
different from the Barker's theory [25]. 



4 Chameleon cosmology 

We begin with the BD chameleon theory in which the scalar field is coupled non-minimally to the matter field via the 
action 

S = J d^x^g-^R-jg^d^dvQ-Ui^+Ji^L,,^, (68) 

where R is the Ricci scalar curvature, is the BD scalar field with a potential U((j>). The chameleon field is non-minimally 
coupled to gravity, ft) is the dimensionless BD parameter. The last term in the action indicates the interaction between the 
matter Lagrangian L m and some arbitrary function 7(0) of the BD scalar field. In the limiting case J(<j>) = 1, we obtain the 
standard BD theory. 



The gravitational field equations derived from the action (68) with respect to the metric is 



(69) 



The Klein-Gordon equation (or the wave equation) for the scalar field is 



1 



(2fl> + 3)D0 = r )+((j)U 4 -2U). 



(70) 



Similarly the energy conservation for the cosmic fluid is 
We shall use the equation of state (EoS) for the fluid p~yp, thus (71) yields p = p (f + to) ~ a . 



(71) 



4.1 Matter collineations 



By defining 



• j 



(72) 



as it is observed it is only necessary to calculate the first component, the rest of them have been already calculated in the above 
sections. 



If Ti =F(0)0-%-, then 
yields 

algebra brings us to get: 



tJ^'pfy -tpj<j>' + tp'j<j> + 2pj<j) = 0, 
W p' ? 2 



J -2 



J P (j) t 

where/' =J ( j ) <p', but we do not obtain more information about the behaviour of some of the functions or J(<j>) ■ 



(73) 
(74) 
(75) 



4.2 Lie Groups 



We need to study the following ODE 



f+ht- x 4' = c[j--4JAt- a +K$v- 



where ' C = "'iW > K =JT2l5> h = const ' 



8a(l-3y) 
3+2<a ■ 

Therefore the standard procedure brings us to outline the following system of PDE 



lit 2 (^,-^) + 277,.-^,-3 



Cr a (j--<j)J^ ) +2U-<j>U$ 



= o, 
= o, 

= o, 



(76) 



(77) 
(78) 

(79) 
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r]„ +hr% +2 Uut -u- cr a (j- l -<$>j^j )§, + acr a - 1 (j - ^j, 
2V - ^ + cr a (j-^yjrif - fa - w„) - 

As we already know, from the MC approach, it must be verified the relationship 



U^-tpU^ 77 =0, 



algebra brings us to get 



Ct 



J(<P) _ t -2+a 



1 



at l ^-2^ t + T] <l) )--fa-(j)J^)T] 

Ui - w„) n + (2u- {-2% t + 77^ 



= 0, 
= 0. 



(80) 



(81) 

(82) 
(83) 



For example the symmetry | = t , 77 = n<f>, brings us to obtain the following restriction on U ((j)) . From Eq. (83) we obtain 
the ODE 

^=K 1_ ^ +2 (H^ (84) 

whose solution has been obtained in section 3, i.e. U (<j)) — Uo<p" <% "~ 2 \ From Eq. (82) we get the next ODE for J : 



'a + 2n-2\ ^ fa+n-2\ J 
Jd>d> = -; 2 



finding that the most general solution is 



\ n J (p 2 

j (<!>)= Jot » {n+a - 2) +c 2 <l> 2 . 



(85) 



(86) 



Theorem 3 The scaling and in particular the self-similar solution admitted for the FE (41-42) have the following form 



<j) = $ (t + t Q )" => U((j)) = Uo^ ( "- 2) = UQ(t + to) 
with p = p Q (t + to)~ a , a = (1 + y)h, while 

/(</>) = /o^ ( ' ,+a - 2 W (f + ro)" +a - 2 . 



n-2 



(87) 



(88) 



5 The Kantowski-Sach model 

We start by considering the Killing vector fields (KVF) 

Z\ = d y , Z2 = cotzcosydy + sinyd,, Z3 = — cotzsinydy + cosyd-, Z$ = d x , (89) 

such that 

[Z 1 ,Z 2 ]=Z 3 , [Z 2 ,Z 3 ]=Zi, [Z 3 ,Z 1 ]=Z 2 , [Z 4 ,Zi}=0, (90) 
in such a way that the metric takes the following form 

ds 2 = -dt 2 + a 2 {t)dx 2 +b 2 (t) (sin 2 zdy 2 + dz 2 ) , (91) 



We find that the metric (91) admits the following HVF 



a 



H = {t+to)d t +( l-(t + t )-\xd x + n-(t+t )j)ydy+ n-(t+t )j)zd z , (92) 



10 



where the scale factors behave as 

b(t)=b (t + t ) and a(t) = a (t + t ) m , Vmet + . (93) 

Tanking into account these result we shall calculate some exact cosmological solutions. It is a straightforward task to carry 
out the calculations for this reason we only show the results. The detailed exposition of the followed method may be found 
for example in [27]. 

5.1 Vacuum solution 

There is no SS vacuum solution for this model 



5.2 Perfect fluid model 

The stress-energy tensor is defined by, Zy = (p + p)u(Uj — pgij, and where we are taking into account the conservation 
principle so, T-j = 0, and the equation of state p = yp , (y '= const ). We find that 

a(t) = a (t + t )^, b(t)=b (t + t ), p=p {t + t Q y 2 , (94) 
and this solution is only valid for 

Y c = l-V2m -0.41421356. (95) 

Note that our solution accelerates since 

q=\-\ = - 1< 0. (96) 

h 2 + V2 

5.3 Time- varying constant model 

In this model we consider the constants G and A as time varying functions in such a way that the FE are: 

G t j = G(t)T tJ - A(t) gij , (G(t)T u - A(t) giJ y j = 0, (97) 

with the constrain T-j = 0. 

We have found the next results 

a(t)=a Q (t + t Q )^, b(t) = b (t + t ), p=Po(f + 'or (r+1)(2W5) 
G = Go(r + f ) (r+1)(2W5 K A = A (f + r r 2 , (98) 

where 

{decreasing Vye [—1,7 C ) ( negative Vy e [— 1,7 C ) 

constant if y = y c Aq « < vanish if y = y c (99) 

increasing Vy 6 (y c , 1] [ positive Vy e (y c ,l] 

where y c is given by Eq. (95). This solution is valid for all EoS y, i.e. there is no restrictions on the y parameter. As above 
q<0. 
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5.4 JBD model with CC 



For this model (this corresponds to the exposed one in section 2) the FE are as follows 

Rij \ 8,iR = +A W*y + - ^<A/f ') + \ (f 

(3+2ffl(^))D^ = 8wr-£D^ iZ ^' , -2^(^-A(^)) J 
where T = T- is the trace of the stress-energy tensor. The gravitational coupling G e s(t) is given by 



(100) 
(101) 



/2o) + 4\ G 



Geff(f) \2co + 3jm 
From the stated theorem of section 2 we already know that the physical quantities behave as follows: 

= o O + *o) n , ==► A(</>)=A <^ ( " +a) =A (f + f r ( " +a) , with « + a = 2. 



(102) 



(o((j)) = const., and p=p (£ + ?o) 



a = (l + y)/i, h = 2 + m. 



We have found the next solution. The scale factors behave as 

a(t) =a (t + t ) m , b(t) = b (t + t ), 

with 



m = ^-{\-y-A)€ [1.2361,4], 



Vye 



-1,1 
'9 



m y= \ = 0, 



withA = \j9y 2 - 10y+ 1. mis not defined Vy e (5,1) • 
where 

« = ^(-i-3y 2 +A(i + r)) e 
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1" 


-T' 2 „ 


, Vye 





«>0, Vye[-l,y c ), 



n 7c = 0, 



« <0, Vye y c , 



1 



«r=i = -2, 



note that y c . is given by Eq. (95). 

l+4y 2 -7y+A(2y-l) 



00 



y(-l+3y+A) 



-3,- 



11 



Vye 



-1,; 



1 



1 



o <OVye [-i,Yi), 0o r ,=O, O >°> Vye (y^- 
where yj = -0.1708203932. Therefore this solutionis only valid if y e (71,5] andy=l. 



-"07=1 



= 1, 



A(<p) = A (t+t y 



A = A (7,0) , <» (0) = cons* = 10 4 , 



(103) 



(104) 



(105) 



(106) 



(107) 



(the recent value of ft) (0) has been obtained from [28]) where the performed numerical analysis shows us that: Ao > Vy e 
(-l,y 2 ), Ao 72 = 0, Aq > OVye (y 2 ,y 3 ), A 0Tc =0, Ao > OVye (y 3 ,i] , and A 0y=1 = 1, where y 2 = -0.4142736105, and 
y c = -0.4142135624. Note that ft e (y 3 , 5] , and then A 0ri > 0. 



p = po(t+t y 



a = (l + y)(m+2), p = p (y,ft>), 



(108) 
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p >OVye (-l,y 4 ), Poy 4 =0, Po < OVye (y 4 ,y 5 ), p 0r5 =0, 

f 



p >OVye (r 5 ,ri), p 0ri =0, p <0, VyG I ft,- 
where y 4 = -0.4187505363, y 5 = -0.4097441700, and y, = -0.1708203932 



p 0r=i = -795.8940815, 



(109) 



Therefore this solution is unphysical, since O > 0, Vy e (yj , |] but p < 0, Vy £ (yj , 5] . Notice that nty c = V2, = 0, 
and Aoy r = as in the above solution. 



5.5 JBD model with potential 



For this model (which corresponds to the particular model, m = 1 , exposed in section 3) the FE are as follows 



g,j = ^ Tij + ^ ( f ^ " + 1 ° 







{2, + 2(o{^))U^ = %%T-(0^^' l + ^U^-2U^), 
where T = 7} 1 is the trace of the stress-energy tensor. The gravitational coupling G e g (?) is given by 



(110) 
(HI) 



Geff(f) 



/2co + 4\ G 



\2a> + 3 J 0(f) 



(112) 



= o (f + r o )" = [/o4>» ( "~ 2) = t/o(f + fo) ( "~ 2) , with « + a = 2. G eff «4r\ 

CO (0 ) = const. and p = p (f + fo)~" , a=(l + /)A, h = 2 + m 



We have found the following solution. The scale factors behave as 

a(t) =a (t + f )"\ b(t) = bo(t + to), 



with 



m = ^-(1-y- A) G [1.2361,4], Vye 
zy 



-a 



, WIy = l = 



with A = ^/9y 2 - 10y+ 1. m is not defined Vy € ( A , 1 ) . 

= </) (r + f )", Geff^r^G^f+^o)"", 

where 

1 . „ r 14 1 

Vye 



3 ' 



n = l(-l-3y 2 +A(l + y)) e 
» >0, Vye [-l,y,), «y r =0, «<0, Vye (y c ,i 



-4 



note that y. is given by Eq. (95), while 

4y(3y- 1 +A) 



5y 2 + 4y-l+A(y+l)' 
Therefore this solution is only valid if y e [— 1 , k] and y = 1 



0o >0, Vye 



-■5 



n r= i = -2, 



</>0y=l = 1> 



£/(0) = t/ o (f + f o ) 



n-2 



{7q = t7o (y, £») , a>(0) = consf = 10 4 , 



(113) 



(114) 



(115) 



(116) 



(117) 



(118) 
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C/ >OVye \(y 2! y e ), t/ 0r =i = 1, 

c/ 0r2 =o, f/o<ovre(r 2 ,r c ), t/ 0y3 =o, 

where y 2 = -0.4142736105, and y c = -0.4142135624. Note that (n - 2) < 0. 

P = Po(t + toy a , a = (l + y)(m + 2), p = p (y,co), 
1' 



p <0Vre -1, 



\(74.r5). Po y=1 



-795.8940815, 



Po r4 =0. Po>OVye (y 4 ,y 5 ), p 0r5 =0, 
where y 4 = -0.4187505363, y 5 = -0.4097441700. 



(119) 




Figure 1: JBD model with potential. Solution Vy € (y 4 , y 5 ). p is plotted in red color. Uq in blue and n in magenta color. 

Therefore this solution is only valid Vy e (y 4 , y 5 ) 9 y £ ., where p > and <j> > while Uq > Vy e (y 4 , y 5 ) \ (y 2 , y c ), 
see fig. (1). Notice that G e ff (y) is decreasing if y € (y 4 , y c ) , constant if y = y c and growing if y € (y c , y 5 ) . As in the above 
solutions we have that that triy = V2, riy c = 0, and Aoy c = 0. We also emphasize that q < 0, Vy € (y 4 , y 5 ) . To end we have 
calculated the values of m in y 4 and y 5 , they are: m n = 1.4117477876367817948 and m rs = 1.4166732719815729598. These 
values will acquire a complete sense in the next model. 



5.6 Chameleon JBD model 

This model corresponds to the exposed one in section 4. In this case FE read: 

J ((b) co ( 1 A 1 / \ U ((b) 

Glj = -f-T tJ + - 2 - -*y W'J + - (f y - ftP*) + -^y, (120) 

(2o) + 3)D</> = T (V- - (<j>U$-2U), (121) 

where the main quantities behave as follows 

(b=(b (t + t ) n , U{<!>) = Uo<!>tM-V = Uo {t + t Q ) (n - 2 \ G eS ^(j)-\ 
CO ((b) = const . and p = p (t + t Q y a , a = (l + y)/i, h = 2 + m, (122) 
J(4>)=J ^ =J Q (t + to) a+n - 2 . 

We have obtained the following solution, n = n(m); 

n = -^, (123) 
m— 1 
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in such a way that n = 0, iff m c = ±1.414213562 = ±y/2, we only consider the positive solution, so, m G R + \ {1} . n < if 

m G [0, 1) U (m c+ ,°°) , and n > if m G (l,m t+ ) . 

The rest of the quantities depend on (m, 7, ft)) , so in order to carry out the numerical analysis it is necessary to fix the value 
of ft) and 7. Setting ft) = 4 • 10 4 , we have studied some equation of state, 7 = 1,1/3,0, — 1/3, i.e. the usual ones. 

For 7= 1 then; <j> = <j> Q (t + fo)" , and then, G e g w = (f + fo)~" , where 

. . , . 1a4 ^ (10002m 4 -2m 3 -60007m 2 + 800 10) , 10/n 

00 = 00 («.y= 1,0J« 10 =- /cn m 4 3 qqqq 2 , , r> ' < 124 > 

(5001m 4 — m J — 9999m- — 3m + 2) 

then 

00=0, mi = 1.414184279, m 2 = 2.000024995, (125) 

and O is not defined when 

m ai = 0.1399429180, m„ 2 = 1.414180724, 
in such a way that <p > if m G (m aii m a2 ) . O < if m G (m ai: m\) , and O > if m G (mi,m2) . t/o behaves as 

(m 2 — 2m + 2) 

f/ o = 0o i - -r 1 , (126) 

2 (m — 1 ) 

therefore Uq has the same roots than O and it is not defined when 

m ai =0.1399429180, m a3 = l, m ai = 1.414180724, 
in such a way that t/o > if m G (m fll , m„ 2 ) \ {m a3 } . t/o < if m G (m fl2 , mi ) , and t/o > if m G (mi , mi) . 

Po = Po { m -Y= 1/3, ft) w 10 4 ,7o = 10) behaves as 

(1 0002;7i 4 - 2»z 3 - 40003/7X 2 - 2m + 40006) 
P ° _ ° 1607T (5001m 4 - m 3 - 9999m 2 - 3m + 2) (m - 1 f ' 

then p = 0, iff 

m , = 1.414184279, m 2 = 2.000024995, m 3 = 1.411747788, m 4 = 1.416673272, 
and it is not defined if 

m ai = 0.1399429180, m« 3 = 1, m„ 2 = 1.414180724, 

finding therefore that p < 0, if m G (m ai ,m a2 ) \ {m„ 3 } , p > 0, Vm G (m$,m4) and p < 0, if m G [m^^mi) . Nevertheless 
a careful analysis shows us that p is not defined when m ai = 1.414180724, note that m (h G (1113,1114) . Thus p > 0, Vm G 
(m^,m a2 ) U (mi^m^) , if m G (m U9 ,mi) then p < 0. 

Note that 

ma, = 0.1399429180 < m„ 3 = 1 < m 3 = 1.411747788 < m„ 2 = 1.414180724 < 

< mi = 1.414184279 < m c+ = 1.414213562 < m 4 = 1.416673272 < m 2 = 2.000024995. 

Therefore this solution is only valid if m G (mj,m a2 ) U (mi^mn) , since in this interval p > 0, O > 0- In fig. (2) we have 
plotted this situation. 

In order to clarify and to compare (with the following results) the results we have write in the following table the roots of 
the constants (n,(j) ,Uo 7 p ) : 
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Figure 2: Chameleon JBD model with 7=1 and m £ (mj,,m4). p is plotted in red color. Uq in blue and n in magenta color. 
O is plotted in green color but it appears under the graph of Uq. 





n 


ri 


n 


r\ 


n 






m c = 1.414213562 








mi = 1.414184279, 






U 




mi = 1.414184279, 






Po 


ot 3 = 1.411747788 


mi = 1.414184279, 




m 4 = 1.416673272 



The cases y = 1/3,0 and —1/3, are quite similar. For example if 7= 1/3, the performed numerical analysis shows us that 
the roots of the constants are as follows (compare with the above table) 





r\ 


ri 


>"3 




n 






m c = 1.414213562 




00 




mi = 1.414197143 


m c = 1.414213562 








m\ = 1.414197143, 


m c = 1.414213562 




Po 


m 3 = 1.411747788 


OTl = 1.414197143, 


m c = 1.414213562 


m 4 = 1.416673272 



and therefore we have the following scenario, see fig. 3 where we have plotted the interval m £ (mj,,rri4). Note the analogies 
with regard to the case 7=1. 




Figure 3: Chameleon JBD model with 7=1/3 and m £ (mj,,rri4). p is plotted in red color. Uq in blue and n in magenta 
color. O is plotted in green color. 

In fig. 4 we have plotted in detail the interval m £ (mi,,m c ). 
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0.0001- 



-0.0001- 




1,41425 1,41431 



Figure 4: Chameleon JBD model with 7=1/3 and m <G (mi,,m c ). p is plotted in red color. Uq in blue and n in magenta 
color. <j> is plotted in green color. 

Therefore this solution is valid for any value of 7, while in the above solution it was only valid for a small interval of y c . 
Remember the solution in the last model, whenm y4 = 1.41 17477876367817948 and m 7s = 1.4166732719815729598. 



5.7 Induced Gravity model 



In this model the FE read 



F tf)G u = Ty + (V;f , - \gitf$*\ + (F giPF (*)) + U tf) gij , 

2U(j>=RF^-2U^ 

where R is the scalar curvature and F (<j>) = (j) 2 /4(0. The gravitational coupling G e g-(?) is given by 



(127) 
(128) 



A 2 



For the metric (91) the scalar curvature is given by: R = 2 ( 2_ -|- 2^ + 2^- + -p 
physical quantities behave as follows: 

<j> = <j> Q (t + t Q y, U((j>)=U (j>^"-^=Uo(t + t f n - l \ G eff «</)- 2 , 
fi)(0) =const, and p = p (f + fo)~" , a = (l+y)h, h = 2 + m. 



(129) 



As we already know, the 



We have found the following solution. The scale factors behave as 

a(t) =a (t + f )'", b(t)=b Q {t + t Q ), 

with m = m(n) : 

m = -n + A, A = \J 2n 2 + n + 2, 
where, as it is observed if n = 0, then m = \f2. <j> Q = O (n, 7) behaves as follows 

2n [3n 2 (1 + 7) +n (1 -2A + ©(1 +A) (1 +7) - (2A + 1)7) +3 (7+ 1] 



00 = 



« 3 (2 + ft) - 7(1 + (0))+n 2 (7(2 + ft) +A) - (O - 1) - 7 (A + 1) - 1 
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note that O = 0, iff 7 = — 1 ,and « 7^ 0. With regard to the constant Uq = Uo(n, 7, (0) we have found the next value 



U = 7° (n (A - n) (CO - 2) + (n 2 + l) (co + 1) + 3) . 
co(n- 1) v y ' ' 




The constant for the energy density, p = p (n, 7, o) , behaves as follows: 



p = ^° (0 O (nco—(l +n) + (n—A) (l — n 2 )) - 2« 3 (© + 1) + 2n 2 (- 1 + 2 (A - n) (1 - to) - co) - 6n) . 




Since the solutions depend on the parameters (n, 7, co) then we need to fix them. For example, if we set For (0 « 4 • 10 4 , then 
the solutions only depend on (n, 7) in such a way that given different values to n then we shall may to study their behaviour. 
For n = — 1, CO « 4 • 10 4 we get (j) > and p > iff 7 < — 1 while f/o < if 7 < — 1. We arrive at the same conclusion if 
n = —2. Note that the cases n = and « = 1 are forbidden. Therefore, the obtained solution is only valid if 7 < — 1 . 

6 Conclusions 

We have studied under the self-similar hypothesis the admitted form of the different unknown functions in several scalar- 
tensor theories. By employing the matter collineation (MC) approach, i.e. calculating the Lie derivate of the effective stress- 
energy tensor with respect to an HVF, and the Lie group method, we have been able to state theorems valid for all the Bianchi 
geometries as well as for a flat FRW metric. We have used both tactics, because with the MC in some of the models studied 
we are only able to obtain relationships between the physical quantities. With this tactic we only obtain self-similar solutions. 
Nevertheless, with the Lie group method we are able to obtain the exact form for each of the physical quantities. Furthermore, 
with this approach we obtain scaling solutions which are more general than the self-similar one. 

In the first of the models studied we arrive to the conclusion that <p « (f + ?o)", and therefore G e g « (f + fo) while 
the dynamical cosmological constant behaves as A » (f + ?o) 2 . In the same way we have deduced that the Brans-Dicke 
parameter cobd (0) must be constant. In the second of the models studied, the generalized scalar-tensor model, the dynamical 
cosmological constant is mimicked by the potential U (<j>) , and therefore we have three unknown functions, F (<j>) , Z (<j>) and 
the potential U (0) . We arrive at a very general result which allows us to outline different scalar-tensor models that admit 
self-similar solutions. Actually this result is in agreement with the fact that we may pass from one model to another through 
conformal transformations. As an example we have emphasized three relevant models, the standard scalar-tensor model, the 
induced gravity model and a specific model which is very similar to the scalar cosmological model where G e ff ~ const. In the 
third of the models studied, following the same procedure, we calculate the admitted form for the unknown functions J(4>) 
and the potential in order to obtain scaling and self-similar solutions. 

Once we have established all these results then we study a particular example, the Kantowski-Sach model. The same 
procedure may be applied to other Bianchi models. We begin by showing that this metric admits and HVF. We explore some 
cosmological models and applying the stated theorems we find exact solutions. We show that there is no vacuum solution. For 
the perfect fluid case, within the general relativity framework, we find that the solution obtained is only valid for a particular 
equation of state, y c , which is not strange in this class of solutions, while the exponent of one of the scale factors is irrational. 
This fact is odd since these constants usually are rational ones. Note that the solution is inflationary since q < without any 
necessity to appeal to a DM component. In the third model, where G and A are considered as time-varying within the general 
relativity framework we have shown that the obtained solution is valid for all value of 7. In this case G may be a growing 
or decreasing function finding that it behaves as a true constant only when 7= y c . The dynamical cosmological constant is 
always a decreasing time function but it may be positive or negative and it vanishes if 7 = y c . Once we know how each physical 
quantity works in the general relativity framework then we explore some solutions for the scalar-tensor models. 

We have not been able to find a solution in the case of a Brans-Dicke model with a dynamical A. This result is quite 
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surprising since following the same procedure we have obtained solutions for several Bianchi models. Nevertheless for the 
BD model with a potential we obtain an exact self-similar solution. The numerical analysis carried out shows us that the 
solution obtained is only valid in a small neighbourhood, £ (y c ) , of y c , the critical value of 7, where the energy density 
and the scalar function are positive. The solution shows analogies with that obtained in the subsection where G and A are 
considered as time-varying within the general relativity framework. For example, G e ff ~ (f + fo) " ma y be a growing or 
decreasing function finding that it behaves as a true constant only when 7 = y c . The dynamical cosmological constant, the 
potential U (0 ) , is always a decreasing time function but it may be positive in the interval of definition or vanish if 7 = y e . The 
model also accelerates since q < 0. Trying to generalize this scenario we also consider a chameleon BD model. We show that 
the obtained solution is valid for all values of 7 instead of only in S (y c ) . In the last model studied, the induced gravity case, 
we find that the solution is only valid if 7 < —1 (a phantom scenario). In the appendix we emphasize the fact that the only 
form compatible with self-similar solutions for the BD parameter (Obd (<P) is (Obd (0) = const. Therefore none of the scaling 
solutions admit a variable (Obd (0) • 



A Study of Eq. (11) 



In this appendix we shall study through the Lie group method the Eq. (11) i.e. 



<j> \2 (0 J t 
Therefore, following the standard procedure, we need to solve the next system of PDE 



W^-^ H = 0, (131) 

-W77 + ^(co w + ^(l-^(/.))77+2r> 2 ^+W(/)77^ + 2 7? w -20^^ = 0, (132) 

f- 2 (^,-^)+2W</)- 1 77,+277^-^„ = 0, (133) 

t~\ + ntt = 0, (134) 

where ff=^_ L 

The symmetry % = t , rj = n<j>, brings us to obtain the following constrain on the function ft) (0) : 



-Wn + 4- ( + (oJ\- ^ ] ) ?? + W$rj, = 0, 



i.e. 



2<a V V (0 



2 ft) / 2(0 \ \ ft) / / \2 ft) 

and therefore we get 

(O^fy + ft) 11 ^~ 

so 

2 

a >0tt> = — --T <O(0) = <»o0 , (00,d€ 

ft) 9 

Therefore if = o f", then ft) (0) = ft>o0 5 = (Oot nS . 
If we substitute this result into Eq. (130) then we get 

n(n-l) = -n 2 (^8-lj-n ^> 8 = 0., 
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this means that 

as we already know. 



(0(<j>) = const. 



References 

[1] P. Bernadis et al. Nature 404, 955 (2000); S. Hanany et al. Astrophys. J. Lett. 545, L5 (2000); A. Balbi et al. Astrophys. 
J. Lett. 545,L1-L4 (2000). 

[2] S. Perlmutter et al. Nature 391, 51 (1998); S. Perlmutter et al. Astrophys. J. 517, 565 (1999); A. Riess et al. Astron. J. 
116, 1009 (1998); P. M. Garnavich et al. Astrophys. J. Lett. 493, L53 (1998). 

[3] C. L. Bennett, et al., Astrophys. J. Suppl. 148, 1 (2003). D. N. Spergel, Astrophys. J. Suppl. 148, 175 (2003). M. Tegmark, 
et al., Phys. Rev. D 69, 103501 (2004). 

[4] A. Benoit, et al., Astron. Astrophys. 399, L25 (2003). A. Benoit, et al., Astron. Astrophys. 399, L19 (2003). 

[5] B. Wang, Y. G. Gong, R. K. Su, Phys. Lett. B 605, 9 (2005). 

[6] C. Wetterich, Nucl. Phys B. 302, 668 (1988). 

[7] B. Ratra and J. Peebles, Phys. Rev D 37, 321 (1988). 

[8] T. Chiba, et al, Phys. Rev. D 62, 02351 1 (2000). 

[9] P. Brax et al, Phys. Rev. D. 81, 103524 (2010). P. Brax et al, Phys. Rev. D. 70, 123518 (2004). 
[10] P. Jordan, Schwerkraft und Weltall, Vieweg (Braunschweig) 1955. 
[11] C. H. Brans , R. H. Dicke, Phys. Rev. 124, 925 (1961). 
[12] PG Bergmann Int. Jour. Theo. Phys. 1, 25 (1968). 
[13] K. Nordtverdt, JR. Ast. Jour. 161, 1059 (1970). 
[14] R. T. Wagoner, Phys. Rev. D 1, 3209 (1970). 

[15] V. Faraoni, Cosmology in Scalar- Tensor Gravity. Kluwer Academic Publishers 2004. 

[16] Y. Fujii and K- Maeda. The scalar- Tensor Theory of gravitation. CUP 2003. 

[17] CM. Will Theory and experiments in gravitational physics. CUP revised edition 1993. 

[18] D. La and P.J. Steinhardt, Phys. Rev. Lett. 62, 376 (1989). 

[19] B. Boisseau et al. gr-qc/0001066. 

[20] S. Lee, JCAP 03 (2011)021. 

[21] A. De Felie et al, arXiv:1005.0868v2. 

[22] B. Boisseau, arXiv:1011.2915v2. 

[23] A.A. Coley. "Dynamical Systems and Cosmology". Kluwer Academic Publishers (2003). 
[24] A. Zee. Phys. Rev. Lett. 42, 417 (1979). 
[25] B.M. Barker. Astro. Jour. 219, 5 (1978). 

20 



[26] C.B. Collins. J. Math. Phys. 18, 2116 (1977). 

[27] J.A. Belinchon. Class. Quantum Grav. 26 175003 (2009). 

[28] B. Bertotti, L. less, and P. Tortora, Nature (London) 425, 374 (2003). 



21 



